Ion-Phosphate Mode in the DNA Low-Frequency Spectra 
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Abstract 

The vibrational dynamics of a DNA molecule with coun- 
terions neutralizing the charged phosphate groups have 
been studied. With the help of elaborated model the con- 
formational vibrations of the DNA double helix with alka- 
line metal ions have been described both qualitatively and 
quantitatively. For the complexes of DNA with counteri- 
ons Na^, if +, Rb^ and Cs+ the normal modes have 
been found, and a mode characterized by the most notable 
ion displacements with respect to the DNA backbone has 
been determined. The frequency of counterion vibrations 
has been established to decrease as the ion mass increases. 
The results of theoretical calculation have been showed to 
be in good agreement with the experimental data of Ra- 
man spectroscopy. 

Key words: DNA, ion, vibration, low-frequency spec- 
tra. 

1 Introduction 

Starting with the early research by Watson and Crick [1] 
it has been a well known fact that a DNA macromolecule 
in the natural form is a salt, a complex of nucleic acid and 
metal cations. The ion environment of the double helix 
plays the dominant role in the helical structure formation. 
The concentration of the counterions in DNA solution de- 
termines the macromolecule helix twisting, bending and 
the recognition of the DNA sites by proteins and drugs 



[2-5]. Therefore, metal cations are of paramount impor- 
tance in the processes of DNA functioning in living sys- 
tems. 

Metal ions interacting with DNA may be divided in 
two groups [4,5]. The first group represents the ions be- 
longing to the diffusion atmosphere of the macromolecule. 
Another group is the ions that are directly bound with 
different structural elements of the double helix. The in- 
teraction of the ions of this group is very specific and de- 
pends on their type. Transition and alkaline-earth metal 
ions are bound principally with nucleic base atoms, and 
the cations of alkaline metals being bound with phosphate 
groups of the helix backbone [5] . 

Under ordinary conditions the thermal fluctuations 
cause vibrations of DNA structural elements and counte- 
rions around their equilibrium positions. The vibrations 
of alkaline metal ions with respect to the DNA backbone 
phosphates (ion-phosphate modes) are likely to be visual 
and vastly intensive in vibrational spectra because of ho- 
mogeneity of backbone phosphate groups and interacted 
cations. Determination of the character of these vibra- 
tions is of great interest for understanding the counterion 
role in the DNA helix mobility and conformational trans- 
formations. 

The vibrations of the double helix structure elements 
may be observed in the low-frequency region because of 
the massiveness of DNA atomic groups and a relatively 
weak interaction between them. In fact, the DNA vi- 
bration spectra show the set of modes in the region lower 



then 250 



[6-17] . The low- frequency spectra of DNA 
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may be conditionally divided in three ranges. The lowest 
range (10 - 30 cm~^) characterizes vibrations of the dou- 
ble helix backbone [16, 19]. The DNA frequencies in this 
range depend on conterion type [13, 15, 16], and humidity 
[9,16]. The middle range (60 - 120 cm~^) is generated 
by vibrations of the hydrogen bonds and intranucleoside 
mobility [16, 19]. A mode strongly depending on counte- 
rion type has been observed in the range within 120 - 250 
cm~^ [17]. In this frequency range an internal vibrational 
mode of the nucleic bases can be also observed [14, 18]. 
Therefore, it is important to estimate a possible range of 
the DNA ion-phosphate modes and their dependence on 
ion type. 

To describe the low-frequency vibrational spectra of 
DNA macromolecule the theoretical approach has been 
proposed [19,20]. However, in this approach the coun- 
terion vibrations have not been taken into consideration. 
The problem of ion effect on the DNA low-frequency modes 
has been studied in works [12,21,22]. The results have 
showed that the DNA low-frequency vibrational spectra 
depend on counterions, but the specific ion mode has not 
been found. 

The purpose of this paper is to describe the low-frequency 
vibrations of DNA with counterions, as well as to deter- 
mine a place of ion-phosphate mode and its dependence 
on counterion type. 

2 Model of DNA with counterions 

In frame of the approach [19, 20] the improved model 
of the DNA vibrations is elaborated [23]. The model 
presents DNA as a double chain of the backbone masses 
Too {P02+'20+C5) with connected pendulum-nucleosides 
TO (sugar -f base) . The counterions are modelled as a masses 
(too) bonded to the backbone. The nucleoside bases of dif- 
ferent chains connected by hydrogen bonds (Fig.l). The 
geterogeneity of the DNA nucleosides is important for H- 
bond stretching vibrations in base pairs and for intranu- 
cleoside mobility [19,20]. In present paper we intend 
to find the place of the ion mode in the haul DNA low- 
frequency spectra, and on this stage the geterogeneity of 
nucleosides and the intranucleoside mobility are not con- 
sidered. 

Let us write the energy of the system: 



The expression for kinetic energy of the n-th monomer 
can be written in the following way: 
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(2) 

where M = mo + to + ma', Yn - the displacement of n-th 
nucleoside and phosphate group mass; / = + maV^ - 
the inertia moment of nucleotide with counterion; I is the 
reduced length of the pendulum-nucleoside, and rg is the 
equilibrium distance between ion and phosphate group; 
On - the angle of displacement of n-th nucleoside from 
equilibrium state O^q] Ig — I sin 9 f„ - the displacement 
of n-th ion from the equilibrium state. 

For small displacements the potential energy of the 
monomer link may be presented as a sum of three terms 
describing the energy of H-bond stretching in the pair, 
the energy of torsion motions of nucleoside around the 
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where Kni, and Uni are the kinetic and the potential en- 
ergy of n-th monomer link of the system, and C/„,n-i is 
the interaction along the chain. 



Figure 1: The four mass model with couterion. (a) The 
monomer links of the DNA. I - the counterions that di- 
rectly bind to the DNA structure elements, II - the diffu- 
sion atmosphere of the macromolecule (b) The displace- 
ments of the atomic groups and ions in monomer link 
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backbone chains, and the energy of ion-phosphate bond: 



(3) 



Here a is the force constant of hydrogen bond stretch- 
ing in pairs, (5„ w ls{6in + 02n) + Yin + i^2n is the base 
pair stretching, /3 - is the force constant of nucleoside vi- 
bration relatively phosphate group, 7 is the constant of 
ion-phosphate vibration. 

In frame of the DNA conformational vibrations model 
(Fig. 1) we are looking for the normal modes of the lat- 
tice of double helix with counterions. To compare our 
calculations with experimental data let as consider a lim- 
ited long-range vibrations of optic type that have to be 
seen in a spectra of real objects. In this case, it is enough 
to consider the homogeneous displacements of the system 
monomer links, where the interaction between monomers 
along the chain (J7„,n-i ) may be waved. 

Using the expressions for the energy 121, and 

let us write the Lagrange equations of motion for this 
system. For the conveniens of the following consideration 
we introduce the variables: 



Y = Yi+Y2, y = Yi- Y2, 
= 0i + O2,V = di-O2, 



(4) 



In variables Q the equations of motion for whole system 
brake up into two systems of connected vibrations: 

mL ■■ m„ •• 

Y + — -0 + = -aaiUe + Y); 
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Here is ap = 2a/M, /3o = P/I, 70 = 7/"^a- 

The normal modes of the system (jSJ and © let us 
find in form: q = gge"*"*. From the system of equations 

the equation for frequencies is obtained: 



(70 -Ix) ) 



(7) 



~ (Po -t^ ) = 0. 

Here po — aoMl'^/I + (3o. The solutions of equation ([7|) 
are the modes wi, uj2, and W3. 

From ((BJ the following equation is obtained: 



a;V + ^'(/3o+7o)+7o/3o = 



(8) 



Here p'^ = ^mjU - 1), 7^ = ^^{mHl/MI - 1), = 
l — ma/M—m^ll/MI. The solutions of © are the modes 

a;4, and ljc,: 



^4,5 



(9) 



Let us analyze the modes 1^4,5. Taking into account 
that for all counterion type ma/M and m^ll/MI are 
rather small (about 0.1), formula ® may be transformed: 



^^4 ~ 70, 



(10) 



iol « /3o. (11) 

Expressions IjlOl) and Hill) show that the modes 074 and 
uj^ depend, respectively, on ion vibration constant, and 
constant of backbone vibration. Analyzing the ratio be- 
tween amplitudes obtained from ©: 
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(12) 



(13) 



(70 - Lu'^)mls ' 

it is clear that on the mode ojn the amplitude of ion vi- 
bration (po) is much lager in comparison with the ampli- 
tudes of phosphate (j/o) and nucleoside (770) motions. So 
the mode 0)4 characterizes mostly by ion displacements. 
Therefore we consider that this mode is the mode of ion- 
phosphate vibrations. 

Let us analyze the equation iQ. As seen, the form 
of equation depends on the value of counterion mass 
{ma). Taking into account that for light ions {Li'^ and 
Na^) the ratio ma/M is very small (about 0.05), the last 
term in equation Q may be waved: 



(70 -w^) 



(ao - w^)(po - t^^) 



Mlt 



(ao - Lo'^m/Mf 



0. 



(14) 

One of the solution of (|14|l is oji w ^/Jo■ The amplitude 
of ion vibration on the mode lui is the largest in compari- 
son with the amplitudes of phosphate (Yq) and nucleoside 
(00 ) motions. So we have found that for light counterions 
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the mode lui may be also considered as a mode of ion vi- 
bration. It means that for light counterions in the DNA 
low-frequency spectra there is a twice degenerated mode 
of ion-phosphate vibrations {lui — lu^ y/jo). 

The analysis for heavy ions and Cs+) chows that 
on the mode wi the amplitude of H-bond stretching is the 
largest. So the structure of the DNA low-frequency spec- 
tra is qualitatively different for light and heavy counteri- 
ons. 

For further analysis of the ion mode let us determine 
the constant of ion vibration. Taking into account that 
counterions form the regular structure around the sugar- 
phosphate backbone of the DNA such system may be con- 
sidered as the ion crystal. According to the classical study, 
the energy of ion interaction in crystal can be presented 
by potential: 

where is a Madelung constant, e is the ion charge, e is 
the dielectric constant, r is the distance between charges, 
J5 and h characterize the repulsion between ions. The 
first term describes the electrostatic attraction and the 
second - Born-Mayer repulsion. Expanding the energy 
(|15|l in the row on small displacements from their equilib- 
rium positions ^ = r — tq, V 7^^/2 the constant of ion 
vibrations is obtained: 

Using the expression for the constant of ion vibration Hlt)|) 
and (|10|l the approximate formula for the mode of ion 
vibration is obtained: 

« J-^^(?-2). (17) 

Analyzing H17|l we can see that the ion-phosphate mode 
strongly depends of the equilibrium distance ro, and ion 
mass nia- Thus we can make the conclusion that the ion 
mode has to decrease with increasing of the ion radius: 

LOLi > tONa > > UJRb > ^Cs- (18) 

The ion-phosphate mode is also depended on the Madelung 
constant Ma ■ This constant is the lowest for dipole {Ma — 
1) and is big enough for the NaCl ion crystal {Ala — 
1.748). In DNA the value of the Madelung constant is 
not known, but we can say that the Madelung constant is 
depended on the ion concentration, and in dilute solution 
it will be another then in the saturated solution. It means 
that the ion-phosphate mode has to be depended of the 
ion force of the solution. 



The constant of ion vibration also depends on the 
other parameters (6, e) that will be considered in the fol- 
lowing section. 

3 Constant of ion vibration 

In this paper we study the conformational vibrations of 
DNA with alkaline metals. So the equilibrium distance 
ro is estimated as a sum of Pauling radiuses of ion and 
oxygen atom of the phosphate group (Table 1). The pa- 
rameter h is taken equaled to values that are known for 
ion crystals LiF, NaF, KF, RbF, CsF [24] (Table 1). 



Table 1. The parameters for ion vibration constant 



Parameter 


Li+ 


Na+ 


K+ 


Rb+ 


Cs+ 


ro{A) 


2.00 


2.35 


2.73 


2.88 


3.01 


b{A) 


0.329 


0.316 


0.327 


0.329 


0.313 


e [25] 


5.0 


6.4 


8.2 


9.0 


10.1 


e [26] 


1.3 


1.5 


1.8 


1.9 


2.1 


Ma {e = const.) 
Ma {e [25]) 
Ma (e [26]) 


1.307 
1.045 
1.050 


1.360 
1.067 
1.066 


1.417 
1.099 
1.089 


1.439 
1.114 
1.099 


1.470 
1.136 
1.116 



Let us estimate the dielectric constant. The value of 
this constant depends on the distance between ions and 
the nature of environment. The dielectric constant de- 
termination for the DNA presents some difficulties [25- 
30]. The value of vacuum dielectric constant (e = 1) is 
not suitable for description of electrostatic interaction in 
DNA because in this case the potential near DNA is too 
high [27]. Very often for solving the Puasson-Boltzman 
equeation the dielectric constant of the DNA-solution sys- 
tem is considered between 2 and 4 [27]. Hingerty with 
coworkers developed an approach where the dielectric con- 
stant presents as a function of a distance [25]: 

e{r) 78-77(r/2.5)2exp(r/2.5)(exp(r/2.5)-l)"2^ (19) 

Lavery with coworkers were improved the function of Hingerty 
[26]: 

e{r) = £oo - ^2^[(sr)2 + 2sr + 2]e-^^ (20) 

were Soo — 78, and s — 0.16. There are also exist other 
distant dependent dielectric functions [27-29], but func- 
tions of Hingrty H19I) and Lavery H20|l are most com- 
mon used. The dielectric constants calculated by formulas 
fHHl and 1201) are shown in Table 1. 

Let us estimate the Madelung constant Ala. This con- 
stant characterizes the interaction of one ion with the 
other charges of the system. In our model the charges 
of the phosphate groups with counterions that neutralize 
them are considered. Also it is taken into account that 
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in the solution there exist negatively charged ions of Cl^ . 
This anions arc situated at the distance about 7 A [4] . As 
usual long-range electrostatic interactions in such systems 
as DNA are catted off on distance about 10 A [31] that is 
approximately equaled to the distance between two neigh- 
bor phosphate groups of one polymer chain. That is why 
only the charges of neighbor ion-phosphate pairs and Cl^ 
anions are considered in Madelung constant calculations. 
The following formula for Madelung constant is obtained: 



M„ = 1 



e{rt)rt <^{rp)rp e{ra)ra 



(21) 



Here e{r) is the dielectric constant, rt is the distance be- 
tween cation and C'l^ anion of the solution, Vp is the dis- 
tance between cation and neighbor phosphate group of 
the same polymer chain, and ra - the distances between 
neighbor cations of the same polymer chain (Fig.l). The 
first term in (|21|l describes the interaction between cation 
and Cl~ anion of the solution, the second term describes 
the interaction between cation and neighbor phosphate 
group, and the third term describes the interaction be- 
tween neighbor cations. The values of the Madelung con- 
stants calculated by formula H21|) are shown in Table 1. 

4 Ion-phosphate mode 

Using the parameters shown in Table 1 and formula 1)16(1 
the values of the constant of ion vibrations (7) is esti- 
mated. The frequencies of ion-phosphate vibrations (0^4) 
calculated by expressions (O are shown in Table 2. 

Table 2. The frequencies of ion-phosphate vibration for 
B-DNA (cm-i) 



Method 


Li+ 


Na+ 


K+ 


Rb+ 


Cs+ 


e = 2 




233 


162 


112 


93 


e = 4 


314 


165 


114 


79 


66 


e [26] 


488 


237 


151 


100 


80 


e [25] 


252 


115 


70 


46 


37 


[17] 


237* 


230* 


150 


110 


95 


[14] 




235 









Note: The asterisk * marks the value determined by us 
from the spectra of work [17] 

Analyzing the results shown in Table 2 we can say that 
the best agreement with the experimental data gives the 
value of e = 2 and Lavery dielectric function [26] . In case 
of ion Li'^ the calculated frequency of vibration differ suf- 
ficiently from the experimental data. The reason of this 
situation may be connected with effects of ion hydrata- 
bility. The ion Li~^ has a very small weight (7 a.u.m.) 
and radius (0.6 A) in compare with molecules of water. 



CO (cm-1) 



200-- 



150-- 



100-- 



50-- 




-\-> 



Na 



Cs" 



Figure 2: The dependence of the low- frequency modes on 
counterion type. • - our results, o - [17], □ - [13]. 



Therefore, the influence of water on the lithium vibrations 
may be very large. Our estimation Li~^ frequency when it 
is in complex with one (wf * = 267cm"^), two (209cto~^), 
and three (180cm~^) molecules of water show that the 
frequency value became sufficiently lower. As seen if one 
consider that lithium ion moves in the solution together 
with water molecules the frequency become very close to 
the experimental value. 

For understanding of the low-frequency spectra de- 
pendence on counterion type the frequencies of all five 
branches of vibrations toi, W2, W3, W4, and are calcu- 
lated (Fig. 2). As we can see all low- frequency modes ob- 
tained with the help of proposed approach decrease with 
increasing of counterion mass. The ion mode {uji, W4) 
is the most sensitive to counterion type. Our calculations 
show the splitting of the ion mode for heavy ions that was 
qualitatively obtained above. The theoretical results for 
the lowest mode (^2,5) and for the ion mode (W4) correlate 
with the experimental data. The counterion dependence 
of the H-bond stretching mode (uj^) is not studied exper- 
imentally. 

In the conclusion we can say that the specific mode 
of ion-phosphate vibration was found. It is in frequency 
range 90 - 250 cm^^ and decrease with increasing of ion 
radius. This mode is very sensitive to the counterion type 
and to the ion force of the solution. As the results pre- 
sented in fig. 2 show all low-frequency modes of the DNA 
macromolecule depend of counterion type. 
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